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The state of art in spin glass field theory is reviewed. We start from an Edwards-Anderson-type 
model in finite dimensions, with finite but long range forces, construct the effective field theory 
that allows one to extract the long wavelength behaviour of the model, and set up an expansion 
scheme (the loop expansion) in the inverse range of the interaction. At the zeroth order we recover 
mean field theory. We evaluate systematic corrections to this around Parisi's replica symmetry 
broken solution. At the level of quadratic fluctuations we derive a set of coupled integral equations 
for the free propagators of the theory and show how they can be solved for short, intermediate 
and extreme long distances. To reveal the physical meaning of these results, we relate the various 
propagator components to overlaps of spin-spin correlation functions inside a single phase space 
valley resp. between different valleys. Next we calculate the first loop corrections to the theory above 
8 dimensions, where we find that it maps back onto mean field theory, with basically temperature 
independent renormalization of the coupling constants, thereby demonstrating that Parisi's mean 
field theory is, at least perturbatively, stable against finite range corrections. In the range between 
six and eight dimensions various physical quantities pick up nontrivial temperature dependences 
which can, however, still be determined exactly. Upon approaching the upper critical dimension 
(d = 6) of the model, scaling which is badly violated in Parisi's mean field theory is gradually 
restored. Below 6 dimensions one should apply renormalization group methods. Unfortunately, 
the structure of RG is not completely understood in spin glass theory. Nevertheless, the first 
corrections in 6 — d to e.g. the exponent of the order parameter can still be calculated, moreover 
exponentiation to this power can be checked at the next order. The theory is, however, plagued by 
infrared divergences due to the presence of zero modes and soft modes. Systematic methods (like 
those developed in the O(n) model) to handle these infrared singularities are not yet available in 
spin glass theory. 



1 Introduction 

Parisi's mean field theory (MFT) [1] is generally accepted as the correct solution of the 
Sherrington-Kirkpatrick (SK) problem [2]. The crucial technical assumption that renders 
the SK model soluble is that every spin interacts with infinitely many neighbours. Such a 
situation could arise e.g. in a system in infinitely high spatial dimensions d — > 00 or in a finite 
dimensional system with an infinitely long ranged interaction. While in some extensions 
of the theory (especially in applications outside physics, like combinatorial optimization 
or neural networks) the assumption of high connectivity is well justified, in real physical 
systems that live in low dimensions and have finite range forces MFT can serve, at best, as 
a zeroth approximation. 
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In order to go beyond the SK model and to approach a more realistic situation, one 
can set up an expansion scheme either in the inverse dimensionality {1/d expansion) or in 
the inverse interaction radius (loop expansion). These will work, however, only if the phase 
space structure of the finite dimensional, finite ranged system is similar to that of the fully 
connected mean field model. 

The physical picture underlying Parisi's solution is of the system having a rugged free- 
energy surface, with many, hierarchically organized equilibrium states in the frozen phase. 
However, a general consensus on whether such a structure can survive in real, finite dimen- 
sional, finite ranged models has never been reached. To the variety of approaches in this 
decade-long debate, ranging from scaling [3,4] and phenomenological renormalization [5] to 
large scale simulations [6-13] and to the first perturbative steps within the 1/d [14] and 
the loop expansions [15,16], respectively, there have recently been added the exact methods 
of mathematical physics [17-21]. With this, a fundamental problem, namely that of the 
definition of an equilibrium state in a random system, has been brought into focus. 

The results to be reviewed in this paper have been obtained under the assumption that 
the "many valley structure" of MFT remains relevant for a finite dimensional, finite ranged 
system. We use the formalism of replica field theory [1], calculate free propagators and set 
out to determine the first loop corrections. All our results will be confined to the Ising spin 
glass just below its freezing temperature. 

The technical difficulties we encounter are considerable, so we can fully compute the 
first corrections to the order parameter and to the excitation spectrum only in high (d > 6) 
dimensions. Near d = 6, the upper critical dimension, we have only partial results. We 
also notice severe infrared (IR) problems here. These are manifestations of the extreme 
sensitivity of the system to changes in distant regions: the unusually strong IR singularities 
in our propagators may well tell, in the language of replica field theory, the same story 
as the mathematical results concerning the chaotic response to changes on the boundary 
[20] . In this context we find it remarkable that the worst IR powers are displayed by those 
propagator components that also describe the chaotic response to infinitesimal variations 
in the control parameters [22,23]. 

If these IR problems turn out to be unsurmountable they will destroy replica symmetry 
breaking (RSB) field theory from within. If, on the other hand, recent efforts by Parisi 
and coworkers [24] and also by ourselves [25] succeed in saving the theory through formally 
exact statements like Ward identities etc., we will find it hard to believe that the theory is 
completely devoid of physical meaning. 

We interpret the goal of this paper in a very restricted sense: apart from citing the 
results by Georges, Mezard and Yedidia [14] in a qualitative manner we shall not be able to 
cover the method of 1/d expansion, nor the important set of papers by Parisi and coworkers 
on the finite size corrections to MFT [26], although both of these approaches might find 
their well deserved place under the title of this paper. What we will try to give an account 
of is the construction of the loop expansion for spin glasses below the freezing temperature. 
Within the space available we will have to be fairly sketchy on this, too. We will not 
provide proofs but will only give some hints as to how the results we display can be derived. 
Concerning the results themselves, especially those for the propagators, we will try to be 
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rather more exhaustive, however, and collect the various propagators, their limiting forms, 
etc. that have appeared over the years scattered in a series of papers [27-30,24]. We will also 
try to demonstrate their use for the calculation of short range corrections in a few cases. 

The plan of the paper is the following. In Sec. 2 we briefly review the elements of the field 
theoretic formulation of the problem. In Sec. 3 we give a sketchy account of Parisi's MFT 
with the purpose of fixing notation. We take the first step beyond MFT in Sec. 4: we spell 
out the set of equations for the free propagators and, for the sake of a first orientation, solve 
them under the (wrong) assumption of replica symmetry, thereby displaying the famous 
de Almeida-Thouless (AT) instability [31], and also the two characteristic "mass scales" 
inherent in the system. In Sec. 5 we abandon the assumption of replica symmetry but confine 
our attention to the "high-momentum" region near the upper cutoff where the equations for 
the propagators can still be solved by a simple iteration. The results obtained in this region 
will be used later when calculating the first loop corrections in high dimensions (d > 8). In 
Sec. 6 we investigate the propagators in the near infrared region, i.e. for momenta around 
the "large mass" (which vanishes like the square root of the reduced temperature near the 
spin glass transition) and much larger than the "small mass" (vanishing like the first power 
of the reduced temperature). We show that in this region the complicated set of (ten) 
integral equations for the propagators reduces to a set of linear algebraic equations, which 
can then be solved by elementary means. We turn to the general study of the formidable set 
of integral equations for the propagators in Sec. 7. Their solution is broken down into two 
steps. First, by exploiting the residual (ultrametric) symmetry of the system, we indicate 
how the problem of the inversion of a general ultrametric matrix can be reduced to the 
inversion of a much simpler object that we call the kernel. Next, in the special case of the 
Hess matrix we show that in the vicinity of the transition temperature it has a very simple 
kernel which can be inverted with relative ease and yields the propagators in closed form. 
The singularities of the propagators determine the excitation spectrum of the system, and 
show that Parisi's solution is marginally stable. From the complicated exact expressions 
for the propagators we extract the limiting forms, valid in the far infrared or extreme long 
wavelength limit (where momenta are comparable to or smaller than the small mass scale) 
in Sec. 8. The propagator components are related to some combinations of correlation 
functions. This relationship is established in Sec. 9 where the physical meaning of some 
of the results obtained so far is also analysed. We step beyond the analysis of Gaussian 
fluctuations in Sec. 10 where we derive the first loop corrections to the order parameter 
and to the excitation spectrum above 8 dimensions. By absorbing the loop corrections into 
the coupling constants we map back the theory onto MFT and thereby show that Parisi's 
solution is, at least in high dimensions, perturbatively stable. However, the renormalized 
quartic coupling blows up as one approaches d = 8 from above. Therefore in the range 
6 < d < 8 one has to rearrange the loop expansion. This is done is Sec. 11. At and below 
d = 6 the loop expansion is bound to break down completely. Although the structure of 
the renormalization group, the usual remedy in such a situation, is not known in spin glass 
theory, we show in Sec. 12 that some information can still be extracted from the logarithmic 
singularities appearing here. In particular, we show that to first order in e = 6—d the critical 
exponent (3 can be computed in agreement with the values of other critical exponents known 
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from above T c [32-34], moreover, this power can be checked to exponentiate at the next order 
in e. Sec. 13 concludes the paper with a brief summary. 



2 The elements of replica field theory 

Our starting point is a standard Edwards-Anderson-like [35] model for N Ising spins in d 
dimensions, with a long but finite-ranged interaction: 



(t) ' v ^ 

In Eq. 1 the external field has been set to zero, the summation is over the pairs 
of lattice sites, z = p d is the number of spins within the interaction radius measured in 
units of the lattice spacing a. f(x) is a smooth positive function which takes the value 1 
for x < 1, and decays to zero sufficiently fast for x > 1, thereby cutting off the interaction 
around r ~ pa. Jij are independent, Gaussian distributed random variables with mean zero 
and variance A 2 . We will be interested in the long distance behaviour of the model (1), for 
which the details of f(x) are to a large extent immaterial, so we can choose / according to 
convenience. 

The long wavelength properties near T c can be extracted by studying the associated 
effective Lagrangean of replica field theory [36] which, for an appropriate choice of /, works 
out to have the form: 

* = -^eN 2 -^)i<i 2 +^ee«s:+ 

a,0 P DVJV a,Pn Pi 

U p p p ( 2 ) 

^ ^ ^ Pl P2 P3 P4 " ' ^ ' 



a,0 Pi 



where the wavevector summations are restricted by "momentum conservation" J2i Pi = 
and by an UV cutoff at p ~ 1/pa. r is the reduced temperature measured relative to the 
mean field value T C MF = A of the critical temperature: 

rpMF _ rp 

t will be assumed small throughout this paper. 

The fields & al3 are symmetric in the replica indices a, (3 = 1,2, ... , n and $ aQ = 0. In 
the truncated Lagrangean (2), first proposed by Parisi [1] and used by several authors since, 
we have kept only that particular quartic term which is responsible for replica symmetry 
breaking on the mean field level. The numerical values of the bare coupling constants w and 
u work out to be 1 in the case of the Ising spin glass. By a slight generalization of the model 
we wish to regard these bare couplings as essentially free (positive) parameters, partly for 
book-keeping purposes, but also because they pick up short range corrections anyhow. 

Now we split the field into an equilibrium and a fluctuating part as 
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with 



Qa0 



xf/ 3 " 









The Lagrangean then splits into four terms: 

£ = £ (0) +£ (l) +£ (2) +£ (3) +£ (4) 

defined as follows: 



£(0) 
£(2) 



= iV 



a,/3 o,/3,7 a,/3 



a,/3 



a</3,7<<5 P 

where G -1 is the inverse of the free propagator: 



£(3) = 



(^2 _ 2r _ 2uq 2 ap){5 Kr 5 Kr + ^ 



(4) 
(5) 

(6) 
(7) 
(8) 



(9) 



>/}V 6 



a,/3,7 Pi 



E E + ? E E *pf« (io) 



a,/3 



Pi 



and finally 



£(4) = 



127V 



Pi P2 P3 P4 



(11) 



Q,/9 Pi 



We shall regard as the bare Lagrangean and 

£(/) = £(l) + £(3) +£ (4) 

as the interaction. 

Statistical averages are calculated with the weight ~ e . For example, the expectation 



value of the fluctuation to first order in £( 7 ) works out as 

J [d*]*^e £ J [d*]^^e £(2) (1 + + . . .) 



< $f > 



J[d*]e^ 



J[d^]e £(2) (!+£« + ...) 



<*f£« > (0) +...(12) 



where < . . . >( ) is the average with the weight e £(2) and < v &p P >(o)= nas been used. 
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,a/3 



By a Wick decomposition Eq. 12 can further be written as: 



a', /J' 



(13) 



The quantity 77 in Eq. 13 is 



%'/3' = 2 



2m 



(14) 



The expectation value of the fluctuation must vanish, so we have to impose the condition 



w + ^E E g «7a(p) + |E g «w(p)w + -" = • ( 15 ) 

P 7^ct,/3 P 



The propagator G, as defined in Eq. 9, depends on the exact value of the order parame- 
ter, so (15) is, in principle, a self-consistent equation for q a p. For a long-ranged interaction 
(z » 1), however, the loop corrections are small, as can be seen by changing the scale of 
the momentum so as to set the UV cutoff at |p| = 1: 



1 

N 



E • 

p\<l/pa 



(2 



J\p\<l/pa Z (27r) d 7|p|<l 



(16) 



If we call <j(p) the propagator G with pa absorbed into the momentum 

(G- 1 (p)) = (p 2 - 2r - 2uql p ){5%8$i + 6%6g) - 

- w{S^qps + 8% qfh + 5%q a5 + 5$q ai ) = (p 2 l + M^U^S (17) 

then the equation of state (15) reads as 

2u 



2Tq af3 + w(q 2 ) af3 + —ql p + 



, 1 1 f jd 
+ z(2nrJ lpl<1 a P 



+ • • • = (18) 



where the corrections are of 0(1/ z 2 ). 

Other physical quantities work out similarly, so by expanding in we generate a 
series in powers of 1 / z with G as the bare propagator. Before going into the analysis of the 
loop corrections, however, we have to solve the problem at zeroth order. 
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3 Mean field 



Mean field theory is recovered in the present framework by letting the range of interaction 
go to infinity, z — > oo. In this limit all fluctuations vanish, = 0, and the free energy of 
the system is just given by the constant in the Lagrangean: 

F = -T\im-C^(q aP ) (19) 

n-»0 n 

evaluated at the solution of 

2u 

2Tq al3 + w(q 2 ) a p + —ql p = . (20) 

Because of the replica limit, n — ► 0, involved in the formalism, one has to solve Eq. 20 
in the space of x matrices which implies that one is able to parametrize the matrix 
q a p in such a way as to permit the limit n — > to be taken. The obvious Ansatz, due 
to Sherrington and Kirkpatrick [2], is to assume that q a p does not depend on the replica 
indices, i.e. 

In the limit n — > Eq. 20 would then become 

2rq - 2wq 2 + ^q 3 = . (20') 

However, the assumption of replica symmetry led to paradoxical results at low temperatures, 
and was subsequently shown to be unstable against replica symmetry breaking fluctuations 
by de Almeida and Thouless [31]. 

After some unsuccessful attempts by various groups, the solution which, at least within 
the context of mean field theory, is now generally accepted was found by Parisi [1]. Parisi's 
solution is based on a hierarchical replica symmetry breaking pattern which reflects the 
ultrametric organization of equilibrium states in the long range spin glass [37]. 

We do not need to go into the details of the Parisi solution here, so we discuss it only 
to the extent necessary for fixing notation. 

We call the sizes of the Parisi blocks p r , r = 1,2, ... R with R, the number of RSB steps, 
going to infinity at the end. For the sake of uniformity of notation it is convenient to add 
Po = n and pr+i = 1 to the two ends of the series p r . The value of the order parameter on 
the r th level of hierarchy will be called q r , r = 1, 2, . . . R. Upon analytic continuation in n 
the series p r becomes monotonically increasing. 

A useful concept we shall make frequent use of in the following is that of the overlap 
between replica indices: a n (3 is essentially the inverse of q a p, that is 

a PI (3 = r if the corresponding q a p = q r . (21) 

Accordingly, the allowed values for a n (3 range from a n (3 = (corresponding to the 
outermost region in Parisi's pattern) to a n (3 = R (in the innermost block). By extension, 
we add a PI f3 = R + 1, corresponding to the diagonal, a = (3. 
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The overlap a D (3 is a kind of hierarchical codistance between replicas a and (3. The 
metric generated by the overlaps is ultrametric: whichever way we choose three replicas a, 
(3, 7, either all three of their overlaps are the same aP\ (3 = aP\~i = (3P\~f, or one (say aC\0) 
is larger than the other two, but then these are equal (an/?>an7 = /?ri7). 

It is evident that any quantity / constructed of the c/'s and depending on two replica 
indices (like {q 2 ) a p = e.g.) depends only on their overlap: f a p = f(a n (3). 

Furthermore, any quantity / Q/ 3 7 depending on three replicas depends only on the three 
overlaps, and since of these at most two can be different, / a( g 7 is, in fact, a function of only 
two variables, e.g. a n (3 and the larger of the other two: 

fafa = /(an/3,max{an 7 ,/3n7}) . (22) 

In the following we will also have to deal with quantities depending on two pairs of 
replicas, like the propagator G a p^s, for example. Ultrametrics implies that of the six 
possible overlaps between a, (3, 7, 5 at most three can be different, which corresponds to 
the elementary geometric fact that a tetrahedron built of equilateral and isosceles faces can 
only have three different edges. For these 4-replica quantities we proposed the following 
parametrization 

— u max{an7,an5},max{/3n7, / 3n<5} ^ Z0 / 

in [28]. This parametrization is redundant in that, according to what has just been said, of 
the four variables displayed in Eq. 23 at least two must coincide, on the other hand it has 
the merit of being symmetric. 

We shall keep to this parametrization throughout most of the paper, in order to ensure 
consistency with previously published material. In Sec. 7, however, the redundancy would 
mask an important property, so there we will make use of the following observation: The 
only case when the two lower indices in Eq. 23 carry independent information is when the 
two upper indices coincide and both lower indices are larger than the common value of the 
upper ones. In all other cases the smaller of the two lower indices is a dummy variable 
which we may drop and keep only max{a n 7, a n S, [3 n 7, (3 n 6} as the sole lower index. 

A frequently used symbol will be 

5 r =p r - p r+ i . (24) 

In the limit R — ► 00, q r goes over into a continuous, monotonically increasing function 
q(x), which turns out to have a breakpoint x±, beyond which it is constant. Although 
for large R the precise choice of the block sizes is largely immaterial, it is convenient to 
arrange the p r 's so that they fill the interval (0, x\) with S r , r = 0, 1, 2, . . . , R — 1, becoming 
infinitesimal, of 0(1/R), and the last one Sr = pn — pr+i = pr — 1 — > (x\ — 1) staying 
finite. 

Returning now to the equation of state (20), under the Parisi parametrization it be- 
comes: 



2rq r + w 



r R 

6t 1t - Pr+lQr + 2 Qr ^ StQt 
t=0 t=r+l 



+ f^ = , (25) 
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or, for R — > oo: 



2rq(x) — w 



f X dtq 2 (t) + xq 2 (x) + 2q(x) C dtq(t) 

JO Jx 



2u o. . 



(26) 



This is easily solved by repeated differentiation yielding the well-known result [1] 

w 



q(x) 



2u X > 



X < X\ 



w 
2u 



(27) 



x\ = qi, x\ < x < 1 



with the breakpoint x\ to be determined from the condition on q\ 

r — wqi + uq 1 = . 



(28) 



To leading order in r, q\ = t/w which could have been read off from Eq. 26 directly, 
by noting that for x\ ~ q\ ~ r all the terms in Eq. 26 are of 0(r 3 ) except the first and the 
one coming from the upper end of the f3 integral. This type of approximation, consisting 
in dropping all replica integrals and keeping only the contribution coming from the vicinity 
of x = 1, will be used frequently in the following and will be referred to as the "innermost 
block approximation". 

The solution given in Eq. 27 is valid as long as r > (T < T c ). With T -> T c , q(x) -»■ 
and one enters the paramagnetic phase. 

The physical meaning of the order parameter q(x) depending on a continuous variable 
had remained a mystery until Parisi showed [38] that the derivative of its inverse 



(29) 



is nothing but the probability distribution of the magnetic overlaps q a b between the equi- 
librium states a and b of the system: 



P(q) = <5(q-q ab )> 
1 \- 

Qab = J^}^ <Si >a< Si >b 



(30) 
(31) 



(In Eq. 30 < . . . > is the thermal average and the overbar is the average over the random 
couplings.) 



4 Instability of the replica symmetric solution and mass scales 

In order to go one step beyond mean field theory, we turn now to the study of the quadratic 
part of the Lagrangean which describes free (Gaussian) fluctuations of the order pa- 
rameter. The main task is to diagonalize the quadratic form in Eq. 8, i.e. to invert (17) and 
obtain the free propagator G. The spectrum of free fluctuations (the singularities of G) will 
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also provide a stability test for the mean field solution found above. Multiplying Eq. 17 by 
G we get the following set of equations for the propagator components: 

{p 2 -It- 2uq 2 1& )G a( 3 n5 -w ^ %8G a (3, lt i ~ w ^ q m G a p^ = S^Sps , 

a<(3, 7 < 5 . (32) 

For a first orientation, let us work out the solution in the replica symmetric case. Then 
G a/ 3 n s wm have only three different components 

G a (3,a(3 = G\ , a / (3 

Gap^a-j = G*2 , a, (3, 7 all different, (33) 
G 01(3,78 = G3 , en,/?, 7, <5 all different, 

satisfying the simple set of equations 

p 2 - 2w<? - Y q2 ) Gl + 4w 9 G 2 = 1 

(p 2 ~ y9 2 ) G 2 - wqd + 3^G 3 = (34) 

p 2 + 6wq — G3 — 4wqG2 = 

where we have used Eq. 20' and put n = 0. The solutions to (34), first written up by Pytte 
and Rudnick [39], 

_ 1 A 4^V \ 



p 2 _ 4« g2 ^ p 2 + 2wq - (p2 + 2u ,g - 4^2)2 ^ 

r 1 f 1 2wg i Aw2q2 \ M) 

2 p 2 -fq 2 \2p 2 + 2wq-fq 2 (p 2 + 2wq-fq 2 ) 2 ) 1 ' 



G 3 



1 4iy 2 g 2 
p 2 - ^q 2 (p 2 + 2wq - ^q 2 ) 2 



display the AT instability [31] in that the first factors have a pole on the positive p 2 axis. 
Note that the instability is due entirely to the quartic interaction: if we had set u = in 
the Lagrangean (2), we would have found an acceptable pole structure in G, at least on the 
level of Gaussian approximation. Also note that there are two "mass-scales" displayed by 
Eq. 35: the unstable mode has a "mass" of order r (the small mass), while the stable mass 
is <D(t x I 2 ) (the "large" mass . 



5 The large momentum behaviour 

When the assumption of replica symmetry is abandoned the solution of Eq. 32 immediately 
becomes very difficult. It is, therefore, important to realize that there are instances when we 
do not need to get seriously involved with the intricacies of the symmetry broken solution. 
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In sufficiently high dimensions the loop corrections to the various physical quantities depend 
upon the short distance behaviour of the propagators only, even in the vicinity of T c , and 
for such large values of the wavevector (around the upper cutoff, i.e. much larger than 
either of the characteristic masses) one should be able to get a solution to Eq. 32 by simply 
expanding G for large p 2 , without having to introduce any particular parametrization for 
the order parameter. 

We have to keep in mind, of course, that on passing from the original model (1) to the 
effective Lagrangean (2) we have dropped all microscopic details, so we can certainly not 
expect to retrieve the precise short distance behaviour from the large-p expansion of Eq. 32. 
The qualitative behaviour will, however, be correct and the results will later allow us to 
illustrate some important points about how the renormalization of the coupling constants 
takes place in high dimensions, and also about the characteristic dimensions of the model. 

In order to implement the program sketched above, we need to decompose (32) so as 
to separate terms containing propagator components of the first, second, and third kind 
(G a /3 t ai3, G a p )a ~ n and G a p^s with a, (3,^,5 all different, respectively). Dividing through by 
{p 2 — 2r — 2ug 2 t/3 ) we also wish to write Eq. 32 in the form of a set of Dyson equations. 
Thus we have: 

G aj 3, a j3 = Ga}, a p + W G^^q^G^^ + w ^ G^^q^G^p (36) 

w^o,/3 u>^a,f3 

G a f3,a~{ = w G^ a ^qf3jG a p :a p + IV Y^ ^07,07^7^0/3,00; + 

w^o,/3,7 

t aj9u>aG a f3,aj'y + wG^ a ^q a ^G a 

(3,13-1 (37) 

u)^a,f3,j 

G a /3, 7 s = wG jS ^ s {qasGap^ + qpsG a p,Pj 

+ qceyG a [)^ a jj 

+ qp-yGafi^s) + 

+w G^qusGap^ + w G { ®l jS q UJ G a p^ 5 (38) 

where 

= p 2_ 2T _ 2uq 2 a/3 ■ ( 39 ) 

These equations are solved iteratively by noticing that for "large p 2 " (i.e. for p 2 » r, q) the 
only term that survives is G^ in Eq. 36, so 

Go/W = ^2 + • • • • ( 40 ) 
This is now substituted into Eq. 37 to yield 

r - wq ^ + (a\ \ 

Lro/3,o7 — —^4 1 , (41J 

and so forth. The results one obtains after the first few iterations are the following: 

G a( 3,af3 = j 2 + - + ^ + ^ {(2t + 2uq 2 a/3 ) 2 + w 2 ([q 2 ) aa + (q 2 )^ - 2q 2 a/3 ) } + 
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+^ {(2r + 2uql p f + 6w 2 t ({q 2 ) aa + (q 2 )^ - 2q 2 ap ) + 
+Auw 2 iLllfl + +^w 2 ql p ({q 2 ) aa + (q 2 )^ - 2q 2 af3 ) + 

w 3 ((q%a + (q 3 )^ - 2q a(3 (q 2 ) a(3 ) } + ■■■ (42) 



G Q7)/ 3 7 = + ^ ^4wTq a/3 +w 2 {q 2 ) a/3 + 2uwq a ^(ql y + q 2 ^ + 

+ ~8 { 12wr2 9a/3 + 6 ™ 2 (^ 2 )a/3 + W 3 (Zqapiq 2 )^ ~ Mafi^ + <?| 7 ) + 



+ (<? 3 )a/3 + q al {q 2 )p~f + ^/3 7 (g 2 )a 7 ) + 12u ^ T 9a/3(9a 7 + <?| 7 ) + 

+2uu> 2 ^(g 2 ) a/ 3(^ 7 + q%) + qlpqajqp-y + J2 I^Qa^qpu^j + 

+ 4u 2 wq a p(q^ + q^ + ql 7 q} 7 )] + ■ ■ ■ (43) 
2w 2 

G a p,y6 = -^-(qaSqp-y + qaryQps) H • (44) 

We will use these expressions in Sec. 10 to calculate short range corrections to the equ- 
ation of state and to the propagators in high (d > 8) dimensions. 

6 The near infrared region 

Now we wish to probe deeper into the structure of the propagator and investigate its 
behaviour in the long, but not extremely long wavelength region. More precisely, what 
we mean is that the momentum may be comparable to the large mass but it remains much 
larger than the small mass scale: 

p 2 > r 2 (45) 

We call this region the near infrared region. 

In order to solve the problem for momenta that are comparable to the large mass 
{p 2 ~ r) we have to construct the equations for the various components of the propagator 
explicitly, by writing out Eq. 32 under the parametrization (23) and letting R go to infinity 
at the end. The set of integral equations obtained this way was first published in [28]; we 
display them here for the sake of completeness. 

(f -2r- 2uq 2 {x)) G x '* + 2w j* dtq(t)G% + 

+2wxq(x)G x '* + 2w dtq(t)G x 'J + 2wq(x) C dtG x f x = 1 , (46) 

Jx Jx 

(p 2 + uq 2 {z) + uq\ - 2uq 2 {x)) G X 'J + w J* dtq{t)G% + 



+w f 1 dtq(t)G x ' x - w qi G x ' x + w f" dtq(t)G x ' x + 

Jx Jx 
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+wq(x) \xG x z % + £ dtG% + zG% + f' dtG^ - G^} + w £ dtq(t)G% + 
+wq(z) { jf 1 d*G? f f - G*?} + ««/(*) {xG^ + £ dfGg} = , x < z < 1, (47) 
(p2 + ng 2 _ ng2( ^ ^ + w £ ^ + ^ + w £ dtq{t)G * ; y + 

•"•'/*//) • f y tiG% + f <i/g; : ;.| • 

• «•</!,•) {•'■g- + £ dtcjj + .rr;; : ; + £ dtG*? - g; : ; - G*f} - 

-«;giG^;« = , x<y<l , (48) 
(V + ng 2 _ uq 2 {y) ^ G *,y + w ^ dtg(t) ( G g + ^ + w J* dtq{t)G ^y + 

+«; /' dtg(t)G^ + u;g(y) { f G x f y + yG* y % + f dtG%+ 

J x K.Jy Jy 

+xG x l % + f\tG x 1 >*-G* : *} + 

+wq(x) ^xG^ + j\tG^ t - Gl^ - wq^ll = , y<x<l , (49) 

\ V 2 G X Z % + wj" dtq(t)G% + wq(y) f dtG% + wq(z) \zG% + £ dtG^ + 

+xG%£ + 2 £dt€%£- 2G%%}=0 , z<x,y<l , (50) 

\p 2 G x : y +w£ dtq{t)G% + wq{y) £ dtG% + wq{z) \zG% + £ dtG^ - G^} + 
+wq{x){xG% + j\tG x z '*-G x z >*} =0 , x<z<y<l , (51) 

(/ + uq \ z) - uq \ y) ) G x % + w jT dttft) (G^;* + Gg) + u; £ d^(t)G*f + 
+^(y) { /' dtG x f y + yG x y % + C dtG% + zG x f z + 

Uy Jy 

+ J* dtG x 4 + xG x f x + 2 £ ,1/g;;;; - 2G^} + 

+^(z) {*G£» + jT dtG*f + 2 £ diG^ - 2G^} = , y < z < x < 1 , (52) 

(p2 + ng 2 (z) _ ng 2 (y) ^ ^ + ^ ^ dtg(t) + G ^ +W J Z dtq{t)C *» + 
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+wq(y) \yG^ x + £ dtG% + £ dt (g% + G' Lr ) + zG%+ 

+ f Z dtG% - G X A + wq(x) \xG*>* + f die;;;; + zG% + /" d*G£f- 

Jy ) I J x Jx 

-G^-G^}+wq(z)y\tG^-G^=0 , x<y<z<l , (53) 
(/ + ng2(z) _ ng2(y) ) ^ + w £ dtq{t) ( G xj + G x^ + w J* dtq{t)G * ; y + 

+«, f dtg(t)G£» + wq(y) { /' dtG£* + yG x y % + dtG*£ + zG%+ 

Jx VJy Jy 

+ £ die;;; + xG%% + £ dtG£? - G£ - ( | + 

+^){^+ £d^ : r-G^} + 

+^(z){£dtG^-G^}=0 , y <x<z<l , (54) 

(V + ng 2 ( Zl ) + uq\z 2 ) - 2uq 2 (x)) G x z[ % + 

+A(z 1 ,z 2 ) + A(z 2 ,zi) = , x<zi,z 2 <l , (55) 

where 

A(Z!, z 2 ) = W f dtq{t)G% x + u, P dtg^G^ + 

JO Jx 

+wq( Z2 ) { £ diG^ - G^} + «,g(x) {*G£& + £ + ziG£* + 

+ £ diG*f - G^ 1 } , x < z! , z 2 < 1 . (56) 

In Eqs. 46—56, the mean field equations (27), (28) have been repeatedly used. 

Note that Eqs. 47, 48 and 49 can formally be obtained from Eqs. 55, 53 and 54, re- 
spectively, (putting z\ = z, z 2 = 1 in Eq. 55, and z = 1 in Eq. 53 and 54) which is why 
we did not write them up explicitly in [28]. The reason for which we display them here is 
twofold. Firstly, it is not immediately obvious why Eq. 55 should contain Eq. 47, etc., since, 
although G*f Z2 is, in general, a continuous function of its variables, this continuity does 
not apply when either of the lower variables takes the value 1. (Whenever this happens we 
have coinciding replica indices, and there is no reason to expect e.g. G Qj a )7< 5, with a, (3,^,6 
all different, to be the same as G a p ta g. In fact, they are different, and there is indeed a 
jump in G%>v Z2 at z 1 = 1.) 

Nevertheless, the equation one derives for G^ 22 does go over into that for G X 'J when 
z\ = z, z 2 = 1. 

The other reason for us to display Eqs. 47, 48 and 49 here explicitly is that the "inner- 
most block approximation" we are just about to apply to Eqs. 46—56 does not commute 
with the limit z — ► 1. 
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Let us now assume that (45) holds. We can than neglect all the 0(r 2 ) terms in (46)- 
(56) which means dropping all the uq 2 terms, but also the terms in xqG or q J G, except for 
the contributions coming from the upper end of the replica integrals (the innermost block 
approximation). For example, the penultimate term in Eq. 46 becomes 



2w f 1 dtq(t)G%$ -► 2wq l G x l ] 



where we have used the fact that (apart from the jump exactly at z±,z 2 = 1 in the lower 
indices) whenever an overlap variable goes beyond the breakpoint x\ it gets stuck to it. 

The simplification gained by dropping the 0(r 2 ) terms is tremendous: the set of integral 
equations (46)- (55) reduces to the following set of linear algebraic equations: 

(p 2 - 2wq 1 )G x 1 'J + 2wq 1 G x l % i + 2wq(x)G x ' X1 = 1 , (46') 
P 2 G X : X + w qi (G** - G x : x ) + wq(x) (G x ^ + G x ' Xl + G x >* x - Gfg) + 

+wq(z) (G X { X X1 - G x { x ) =0, 1 > z > x , (47') 
p 2 Gl% + w Ql ((;:■<>,. - (;[':) + wq(y) (g x % + CI:',?) + wq(x) (G x y x xi - G X 'J + 

+ G X > X X1 - G x ; x ) = 0, x<y<l , (48') 

p 2 gt v x + mi ((■:■",■ - <%%) + m(x) ((;■:■■!.. - (•::■[) + wiv) (g x : xi + (■[";: ■ 

+ G x 1 ' x xl -G x i: x 1 )=0, y<x<l , (49') 
p 2 Gfy + 2wq(y)G*? + 4wq(z) (c%* x - G x f x ) =0 z < x, y < 1 , (50') 
p 2 G x % + 2wq(y)G^ + 2wq(z) (g x ^ x - Gfg) + 

+2wq(x)(G*f xl -G x z: x )=0, x<z<y<l , (51') 
p 2 G x : y + wq(y) (Gffi + + 2G x x[ * x - 2G x f x ) + 

+2wq(z)(G x >y x -G x : £)=0, y<z<x<l , (52') 
p 2 G x % + wq(y) ((;■:■;:. + G x ' Xl + G%f x - Gfg) + wq{x) + G x %- 

-G x : x -G x : x )+wq(z)(G x * x -G x 1 f)=0, x<y<z<l , (53') 
P 2 G X % + wq(y) ((I'-;; + G x ' Xl + G x >* x - G x f x + G£* - G x z > x ) + 

+w(q(x) + q(z))(G x % 1 -G x x *)=0, y<x<z<l , (54') 
p 2 G x z[ % + wq(z 1 ) (G X > X X1 - G% X A ) + wq{z 2 ) (G x z[ x Xl - G£*) + 
• "•'/!•'•) (G x z[ x x + G%% - Glf + G°g + 

+ G x x[ %-G x ^)=0, x<z 1 ,z 2 <l . (55') 
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In Eqs. 46'— 55' we have set r = wq\, which is the correct stationary condition at this 
order. 

It is evident from Eqs. 46'— 55' that in the given approximation each component of G 

.1 wqi wq(x) wq{y) wq{z) 
is times a function of — 5-, » — , k — and k — : 

pZ pZ pZ pZ pZ 

G=- 2 g{^) , P 2 »r 2 (57) 

pZ pZ 

where q stands for q±, q(x), q(y), q(z), according to the component in question. 

The solution of Eqs. 46'— 55' is particularly simple if all overlap variables are greater 
than x\. The propagators in this limit will be seen to describe fluctuations in a single phase 
space valley. Then the set (46')— (55') collapses to 3 simple equations with the solutions: 

n x llXl _ 1 (. , Iwqi Aw 2 q\ \ 

G i,i ~ 3 1 + ~2 , o...„ + 




p 2 + 2wqi (p 2 + 2wqi) 2 I 

r n - n - I 1 2wqi 1 4w2g i ^ (^) 

1>X1 - 2 ' 2 p 2 + 2w qi (p 2 + 2w Ql ) 2 J 1 j 



4w 2 g 2 



pi (p2 _|_ 2y;g 1 ) 2 



Comparing these with the replica symmetric propagators given in Eq. 35, we see that 
to the order regarded here (i.e. neglecting r 2 terms ), they are the same. The difference is, 
of course, that Eq. 58 is supposed to be valid only for p 2 3> r 2 , so for the time being we 
keep away from the region where the instability showed up. 

Two particularly illuminating combinations are the transverse and the longitudinal 
propagators 

r — r iXl ' X1 oriu^i 1 «i,n _ ^ 

u ± — z,u l,x 1 ~ ru x 1 ,x 1 — ~2 > 

P ~ (59) 



G\\ = 0^-40*^+30^ = ^^ 



which are precisely of the form of the Gaussian propagators in a massless phase. 

Now we display the complete set of propagators valid for generic overlaps in the near 
infrared limit p 2 3> r 2 : 



2wqi 2w 2 q 2 2w 2 q 2 (x) p 4 + 8wq\p 2 + 8w 2 q 2 
rj J ' ' ' ' ' ~* 7 {p 2 + 2w qi ) 2 



8w q (x) 



p A {p 2 + 2wq\) 2 



(60) 



G x,x = 1J( 1H 4 ^i£ l w 2 q 2 w 2 q 2 {x) 8w 3 q 2 {x)q(z) 

M p 2 } V p 2 J p 2 (p 2 + 2wq\) 2 p A p A p 2 + 2wq\ 




G\ 



X 



X 



rix,y 

U 2,2 



r*x,y 
^ z,x 



n%,y 

^ z.x 



r<x,y 
w z.x 



nx,y 

^z,z 



S~1X,X 
^Zl,Z 2 



+ 



Aw A 2q\x)+q 2 {x)q 2 {z) 
p 4 (p 2 + 2wq 1 ) 2 



X < z < 1 



r y,x = wq(x) \ 2wqi _ wqxp 1 + w 2 q{ wq{y) _ Aw 2 q 2 {x) +q 2 (y) 
" 4 1 + " 2 (p 2 + 2w qi ) 2 p 2 p 2 p 2 + 2w qi 



P 1 



V 



Aw 3 q 3 (y) + 2q(y)q 2 (x) 



+ 



p2 ^2 _|_ 2wqi) 2 
wq(y) 



x < y < 1 



2wqi wqip 2 + w 2 q 2 wq{x) \w 2 q 2 (y) + q 2 (x) 



P* 



(p 2 + 2wq\) 2 p 2 p 2 p 2 + 2wqi 



+ 



Aw 6 q (x) + 2q{x)q 2 (y) 
p2 (p2 _|_ 2wqi) 2 



< X < 1 



w 2 q 2 {z) (p 2 + 2wq\ — 2wq{x)){p 2 + 2wq± — 2wq{y)) 



{p 2 + 2w Ql ) 2 



^l(p 2 + 2 Wqi -2wq(y)) 



z < x, y < 1 
wq(z) w 2 q 2 (x) + w 2 q 2 (z) 



p" 



p 2 + 2wqi (p 2 + 2wq±) 2 I ' 

x < z < y < 1 



2wq(x) 



W q(y) + W q( Z ) - 2 W 2 g!M±M^)±^M + 



+2w' 



2wq{y) 



3 q 3 (y) + 2q 2 (x)q{y) + q{y)q 2 {z) 



(p 2 + 2wqi) 2 



p 2 + 2wqi 

x < y < z < 1 



p" 



+2w 



rty,x 
^z,y 



W q(x) + W q(z) - 2w 2<M±^pM±l^ + 

p 2 + 2wgi 

y < x < z < 1 



3 g 3 (x) + 2g 2 (|/)g(x) + g(a;)g 2 (z) 
(p 2 + 2u;g 1 ) 2 

^(^ + 2^-2^)) 



p 2 + 2u>qi 



w q{y) + w 2 q 2 (z) 
(p 2 + 2w qi ) 2 



+ 



w 2 q 2 (x) + w z q(z!)q(z 2 ) - 



, y < z < x < 1 

Aw 2 q 2 {x)(q(z 1 )+q(z 2 )) 



p 2 + 2wqi 



2w 4 q 2 (x)(q 2 (z 1 ) + 2q 2 (x) + q 2 (z 2 )) 
{p 2 + 2w qi ) 2 



, X < Z\ , Z 2 < 1 



These formulae can be verified by a direct substitution into Eqs. 46'— 55'. 
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Note the simple excitation spectrum displayed by Eqs. 58—69: the propagator compo- 
nents have two poles, one at p 2 = —2wq\, the other at p 2 = 0. When we work out the 
complete expressions for the propagators below, we shall see that these singularities form, 
in fact, two continuous bands which, having a bandwidth of the order of the small mass, 
cannot, however, be resolved in the near infrared. In addition to these cuts, we shall also 
find a pole at p 2 = 0. 

The above formulae, derived under the only assumption of p 2 3> r 2 should also cover 
the region near the cutoff. Indeed, the large p expansion of Eqs. 60—69 reproduces Eqs. 42, 
43 and 44 provided the coefficients of 1/p 2 , 1/p 4 , etc. there are evaluated in the Parisi 
parametrization to leading order in r. Eqs. 60—69 will be used when calculating loop 
corrections for 6 < d < 8. 

7 The complete solution for the propagators 

The expressions derived in the previous Section are valid in the range where the momentum 
is much larger than the smaller mass scale involved in the problem. If we tried to continue 
these formulae into the exceedingly long wavelength limit, i.e. to momenta around the 
small mass scale, p 2 ~ r 2 , or even much smaller, p 2 <C r 2 , we would encounter truly 
unmanageable infrared singularities, with the propagator components blowing up like p~ 6 . 
In order to uncover the true behaviour in this extremely long wavelength limit, which we 
will call the far infrared region, we have to return to the original set of Eqs. 46—56 for the 
propagators. Unfortunately, we have not been able to devise any approximation scheme 
(analogous to the "large p expansion" or the "innermost block approximation" working 
near the upper cutoff and in the near infrared, respectively) that would directly give the 
propagators for the far infrared. Therefore we have to face the task of solving the complete 
set of Eqs. 46-56. 

This task may, at first sight, appear utterly hopeless. Nevertheless, it turned out to 
be possible to find a solution to the set (46) — (56) in closed form (even in a more general 
setting, with an external magnetic field acting on the system) which was published by two 
of us in [28], with a minor error corrected in [40]. 

In retrospect, it is clear that the solvability of the problem depends upon the combina- 
tion of two independent factors. One of them is the residual symmetry which allows one 
to reduce the inversion of not only the Hessian, but also of any ultrametric matrix to the 
inversion of a much simpler object we called the kernel in [40,41]. The other is that close 
to the transition point this kernel turns out to be very simple. These two factors deserve 
a separate discussion in their own right which is what we present in the two subsections 
below. 

7.1 The inversion problem of a generic ultrametrix matrix 

For the sake of the discussion in this subsection we shall assume that both the replica 
number n = po and the block sizes pi,P2, ■ ■ ■ ,Pr, featuring in the Parisi construction are 
positive integers such that Pi+i is a divisor of pi, i = 0, 1, 2, . . . , R — 1. This will allow us to 
stay within the limits of well established mathematics throughout this subsection. When 
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coming to the application of the results in the next subsection, we shall have to consider, 
of course, the replica limit n — > 0, along with the p^s to be continued into the interval [0, 1] 
and R — > oo. 

The free energy F in the MF approximation is a functional of the order parameters q a p, 
with q a/ 3 = qp a and q aa = 0, a, (3 = 1, 2, . . . , n. Since none of the n replicas is distinguished, 
F must depend only on permutation invariant combinations of the q a ps, such as those in 
Eq. 6, for example. This means that F is invariant under the action of the elements II of 
the permutation group S n : F = F(q aj3 ) = F (q n(ahn(/3) ) . 

The physical value of the free energy is obtained from the functional F by evaluating 

OF 

it at a saddle point where — = 0. The solution of the saddle point equations is a point 

OQa,p 

in the \n(n — 1) dimensional space spanned by the independent components of q a p, thus 
the order parameter q a p, which is usually spoken of as a matrix, is, in fact, a vector. We 
do not want to introduce a separate notation for this vector, but propose to think of it as 
a column vector, whose components are the elements (above the diagonal) of q a p listed in 
some fixed order. 

The symmetry of F under permutations does not imply that all the solutions of the 
saddle point equations must respect this symmetry. In fact, we know from de Almeida and 
Thoulcss [31] that at low temperatures the replica symmetric saddle point becomes unstable 
against fluctuations that break permutation symmetry. Symmetry breaking gives rise, like 
at any other phase transition, to the reduction of the original symmetry group to one of 
its subgroups. Parisi's Ansatz [1] is, in fact, nothing but a concrete choice for this residual 
symmetry, one that, as it turned out later [37], embodies the physical assumption of the 
existence of many equilibrium states with an ultrametric organization. For this reason we 
proposed the name ultrametric group for this particular residual symmetry in [41]. 

The vector pointing to the Parisi saddle point is invariant under the action of the 
ultrametric group. The set of \n(n — 1) x \n{n — 1) matrices that, acting on the Parisi q a p, 
perform the permutations belonging to the ultrametric group (and thus do not change the 
form of q a p) constitute an (in general, reducible) representation of the group. 

Now consider the ^n(n — 1) x \n(n — 1) dimensional matrices that commute with all 
the matrices belonging to the previous set, i.e. with the matrices representing the elements 
of the ultrametric group. These are the matrices that we called ultrametric matrices in 
[41] where we also gave a detailed description of their structure. If we transform our 
coordinate system in replica space such as to decompose the representation of the group 
into its irreducible components, we find that, by virtue of the Wigner- Eckart theorem, in 
this new representation all ultrametric matrices will be block diagonal. 

The Hessian or stability matrix (essentially the inverse propagator) is an important ex- 
ample of an ultrametric matrix. Constructing the irreducible representations of the residual 
symmetry group one can therefore make significant progress towards the diagonalization and 
inversion of the Hessian or of any other ultrametric matrix. 

The approach indicated above is the standard one. An alternative that may involve 
seemingly ad hoc steps but leads to perfectly identical results is to directly construct the 
orthogonal subspaces which are closed under the action of any ultrametric matrix and obtain 
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the block diagonal form by transforming the Hessian to this new basis. This is what we 
did in [42], exploiting the consequences of this block diagonalization fully in [41]. (The 
decomposition of the reducible representation mentioned above into irreducible parts has 
since been constructed by Bantay and Zala [43].) 

A most detailed account of our approach can be found in [41]. The same results 
were somewhat later reproduced via purely algebraic means by De Dominicis, Carlucci 
and Temesvari [44] by using an efficient tool called the replica Fourier transform (RFT). 

Whatever the procedure, the central result one obtains is the following. In the new 
basis the Hessian breaks up into a string of (R + 1) x (R + 1) blocks followed by a string 
of 1 x 1 "blocks" along the diagonal. The (R + 1) x (R + 1) blocks are called, for reasons 
of little interest, longitudinal-anomalous (LA), and they are labelled by an index k = 
0, 1, 2, . . . , R + 1. The sector where the lxl blocks appear, i.e. where the transformation 
actually diagonalizes the Hessian, is called the replicon (R) sector. 

Let us call our generic ultrametric matrix (in the original coordinate system) M and 
its inverse G. Their components can be parametrized as in Eq. 23. Let us further call 
the corresponding matrices in the new representation M and G, respectively. As shown in 
[41,42], the diagonal elements of M in the replicon sector are given by: 

R+l R+l 

K = E*»E * W " M t-\s - Mg-i + ACi,.-i) > (70) 

s=k t=l 

where the discrete indices r, k, I needed to label these matrix elements take on the values 
r = 0,1, ... ,R, and k,l = r + 1, r + 2, . . . , R + 1, respectively. Similar formulae hold, of 
course, for Q. 

gZ = X> X>(^;I - - giu + gilVi) • ( 71 ) 

s=k t=l 

The combinations appearing in Eqs. 70 and 71 can alternatively be regarded as (double) 
RFT's [44]; in that context, to conform to earlier notation, it is useful to denote them as 

(72) 

GZ = F Z , r + l<M<l?+l. 

The idea of a transform sharing the convolution property with ordinary Fourier trans- 
forms was first proposed, for the continuum limit, by Mezard and Parisi [45] in the context 
of random manifolds. It was later recognized by Parisi and Sourlas [46] as a Fourier trans- 
form within the p-adic number approach (again limited to R — > oo, n — > 0). In this context 
Parisi and Sourlas could derive, for the replicon sector, the relationship • = 1 which 
in our geometric (or group theoretic) approach follows from the fact that M • G = 1 and 
that these matrices are diagonal in that sector. The extension of the RFT to the discrete 
case and to the LA sector through the use of generalized Parisi boxes pi' ,s \ as described 
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below, was made in [44]. Finally, the connection to standard Fourier transforms over a 
group was given in [47]. 

Turning to the longitudinal-anomalous sector, there we have the (R + 1) x (R + 1) 
diagonal blocks labelled by the index k = 0, 1, . . . , R + 1. The matrix elements work out to 
bef| 

NQ" = k k {r)5f; + \k?5?-V , r, s = 0,l,...,R , (73) 



where Afc(r) is shorthand for 



MT'" , k > r + 1 



I K+\,r+X , k<r + l 

S (l) = p (l) _ p (l) +i t 1 = 0,1,... ,R + 1 , s = 0,l,...,R (75) 

and 

P V = {l s ' S ~\ . (76) 

\ 2p s , 8>l V 7 

We shall refer to the objects and F as the kernel or RFT of M and G, respectively. 
K^' s is given in terms of the original matrix elements as 

R+l 

K k S = Pt' S \ M t' S ~ M t-i) ■ (77) 

t=k 

It is here that we make use of the observation made below Eq. 23, and keep only the 
larger of the two lower indices of M. Formula (77) could not have been cast into the simple 
form above had we kept to the original, redundant parametrization. 

The weight p^'^ is defined as 

p J r,s) = I 2p t r <t<s . (78) 



The corresponding formulas for GV are 




^ = X 7^e + ^ fe - 1) , (79) 



with the kernel of G given by 



R+l 

F r k S =Hpt ,S \ Gr t S - G l-i) ■ (80) 

t=k 



"In order to make contact with earlier notation, we note that the quantity M^'\ we use here was called 
X(r;k,l) in [41]; M^ a here is related to Af',, in [41] by a similarity transformation, and the relationship 
between the quantity Al^ in [41] and Sl k ~ 1 ' 1 here is A?' = ^Si h ~ 1 \ The kernel is, however, invariant, so 
it is exactly the same as the quantity called Kk(r,s) in [41], and the Dyson equation (81) also remains the 
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Now the equation G_ • M = 1 expressed in terms of the diagonal blocks reads 

XX/ x^y + ^'*^- 1) )(e A *(*) + = C - 

or 

1 1 1 fl t <5 (fc_1) t 

Fr = "A^^A^)-A^)g F ^^^ • (81) 

If we regard M as the self-energy matrix and G the propagator then Eq. 81 is just 
the Dyson equation connecting their respective kernels or replica Fourier transforms. This 
Dyson equation was first obtained in [40] for the continuum limit and in [41] for the discrete 
case. 

Given the matrix M to be inverted, K r ^ s can be computed using Eq. 77. Eq. 81 is then 
a set of matrix equations in the unknown Fjj s (one for each value of k). Suppose we are 
able to solve it. Next we have to invert the relation (80) to get G itself in the LA sector, 
and invert Eq. 71 to get it in the R sector. One of the great merits of the RFT approach is 
that these inversion formulae work out to be fairly transparent. The inverse of Eq. 80 is 

aGT'' = E-^(^''-%) ■ (82) 
fc=o P k ' 

that of the double transform in Eq. 71 is 

u 1 v 1 

rG^v= — 2 - - F kJ+i + F k+i,i+i) ■ ( 83 ) 

k=r+l Pk l=r+l Pl 

As for aGu T v we find 

A^ u ,v — A<~r u + A^ v ~ AW i \° 6 ) 

whence the full propagator is obtained as 

or = A ar , 

(84) 

G r ' r = P G r ' r + A G r ' r 

Now we proceed to apply the above formalism to the special case of inverting (17) in 
order to determine the free propagator. 

7.2 The exact solutions for the Gaussian propagators near the transition temperature 

Let us now take the inverse bare propagator for the ultrametric matrix M of the previous 
subsection. As seen from Eq. 17, its components are 

P 2 + M R+i,R+i = P 2 -2r-2uq 2 r , r = 0,1,..., R 

M r ^ +1 = -wq(mm{r,s}), r,s = 0,1,..., R (85) 

M R+i,t = ~ w Qt, r<t = 1,2,..., R . 
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In the truncated model (2) all other components are zero. 

Now we substitute these into Eq. 70 and take the continuum limit, i.e. perform the 
analytic continuation in the p's into the interval (0,1) and let R — * oo, according to Parisi's 
prescription [1]. In the new representation we get for the diagonal components in the 
replicon sector 



= p 2 +uq 2 {s) +uq 2 {t) -2uq 2 (x), < x < s,t < 1, (86) 

where we have also used the equation of state (26). Since in the replicon sector M is 
diagonal, the A^ ) appearing in Eq. 86 are nothing but the eigenvalues of the Hessian matrix 
M^ ) defined in Eq. 17. These eigenvalues fill a continuous band with lower band edge 

2u 

\(°\x;x,x) = and upper band edge A^(0; 1, 1) = 2uq\ = —^t 2 + ■ ■ ■ 

w z 

0<A(°)(x;M)<^r 2 . (87) 

The LA kernel corresponding to the simple ultrametric matrix (85) is obtained from 
Eq. 77 as 



K% = -4wq(mm({s, t}), < s, t < 1 . 

Note that this is independent of the lower variable k. (If we had also kept the Tr(i^ 4 ) 
invariant in Eq. 2, as in [40], we would have a ^-dependent kernel.) 
In the continuum limit Eq. 81 then becomes 



x>y _ 4wq(mm{x,y}) 2w f 1 Ak) , < ,v„ 
Af(x)Af(y) A^(y)Jo 



x,t 
k ' 



< k, x,y < 1, 



where 



A { t k) = (h(k -t) + ^5(k -t) + 9{t - k)j dt (90) 
is the continuum limit of — ^5[ k ^ and|| 

Ai 0, w = {^ + A<ol(I;I ' fc) ; HI PD 

The kernel in Eq. 88 has a remarkable property, namely that it depends only on one 
(here, the smaller) of its upper variables. This is, in fact, the other crucial factor (in addition 
to the symmetry analysed in the previous subsection) that allows one to obtain the bare 

To avoid confusion we note that Eq. 90 is also the continuum limit of the negative of a quantity which 
was, unfortunately, denoted by the same symbol A^ fe) in [41]. 
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propagator in closed form. Indeed, for a kernel depending only on the smaller of its upper 
variables the solution to Eq. 89 becomes a product of two factors, one depending on x, the 
other on y. 

We write it in the following form 



? x, y 4 (min{a;, y}) 3> fc (max{x, y}) 



k ~ Wk A< 0) (min{x, y}) Af (max{x, y}) 



where ^ are two independent solutions of 



d 1 d ^±^_ 



■ ■ , *t( x ) = mr- (93) 

dxwq{x)dx ky 2Aj 0) (x) 

and the Wronskian 

wq(x) 

is independent of x. (The dot means derivative w.r.t. the argument). 

For the details of the derivation of we refer the reader to [28,40], here we content 
ourselves with simply giving the results for Fj°' y . In order to simplify the notation slightly 
we introduce 



-2 — U 2 x l 2 /ncA 

ur Iwqi 



With Eqs. 86 and 91 Af\x) now becomes 

a, w w = ( " 2 i 1 + • x - k (*) 
y p 2 , x > k. 

F^ ,y will be expressed in terms of the solutions Q\ and Q2 of the Gegenbauer equation 

(1 - e)Q = 2G (97) 

belonging to the initial conditions 

Qi = C(0 , C(0) = 1 , C(0) = . . 

g 2 = s(o , 5(0) = , 5(0) = 1. 1 ' 

The variable £ is related to the original overlap variable x by 



x I k 
2p \ Ap- 



1 

2 \ ~2 



t = ~ I 1 + TT2) > (") 



and we will frequently use abbreviations like C(£) = C x or C k = C x= fc etc. . 
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We now have 



2 T^ x >v 
P F k 



2 
P 



1 + k 2 - x 2 l+ k 2 -y 2 A k a k + B klk S k 



Ap 



4p2 



{A k a k + B k -f k )C y 
2 



2 \l 1 + TI? ( a k + ^pk ) C k + 7 fe C fc ] S y 



4p 2 



p 2 F *,v = 



P F k = 



1 



Sr. 



P 1+ k 2 -x 2 A k a k + B k j k S k 



7s/ 



Ap 2 



2 
P 



A k cosh 



x — k 
P 



k sinh ( 



x — k 



\ P 



A k a k + B k ^ k 



where the notations 



A k = 

B k = 

lx = 

<Jx = 



k 2 

2 f + W Sk 



s t + 

cosh 
sinh 





k 2 


4p2 


f Xl — X 




V p 


f Xl — X 




V P 



s k , 

1 — Xl 
P 

1 — Xl 

p 



sinh 



cosh 



Xl — x 
P 

Xl~ X 

p 



,x < y < k < xi, (100) 
x < k <y < xi, (101) 

k<x<y<xi, (102) 

(103) 

(104) 
(105) 
(106) 



have been introduced. 

F^ ,y is a continuous function of x, y and k. For y < x, Eqs. 100—102 apply with x and 
y interchanged. Whenever either x, y or k goes beyond x\, F^ ,y becomes a constant in that 
variable. 

In the above notation the Wronskian in Eq. 94 works out to be 



W k = A k a k + B k -/ k . 
It can be shown that the roots of 

W k (p 2 = -X) = 



(107) 
(108) 



give the LA eigenvalues of the Hessian M*- ^ . 

Writing up Eq. 107 explicitly with the help of the definitions (103)- (106) we can show 
that Eq. 108 has infinitely many positive roots for any given < k < x\. These can be 
labelled by a discrete index as X m (k), m = 0, 1, 2 . . .. The eigenvalues belonging to m = 
are of the order r and, with k varying between and x±, they form a continuous band: 



ar (l-^ T + ...)<A (fc)<2r(l + ^ s r + ...). 



(109) 
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All the other eigenvalues are of C(r 2 ). For a given k they fall off like ~ while for a given 
m they form a continuous band of width of 0(t 2 ) as k is varied. The bands corresponding 
to m = 1,2,... partially overlap and even the largest (Ai(xi)) of them is smaller than the 
upper edge of the replicon band given in Eq. 87. Therefore the small, 0(t 2 ) eigenvalues 
can be regarded as forming a single continuous band spanning the range [0, ^t 2 }. 

Replica symmetry breaking resolves part of the degeneracy of the AT eigenvalues: the 
large eigenvalue is split and smeared out into the narrow band (109), while the small (nega- 
tive) eigenvalue, in addition to being split, also gets shifted, so that the smallest eigenvalue 
is now zero: RSB has cured the AT instability. 

It is clear that these continuous bands of eigenvalues will show up as narrow branch cuts 
in the propagators. In the near infrared limit, with squared momenta much larger than the 
bandwidths, these cuts appeared as simple poles. We have now learnt the precise structure 
of the singularities so we are prepared to go into the far infrared limit. Before turning to 
this we write up the various propagator components with the help of the (continuous forms) 
of the inversion formulas given in Eqs. 82, 83 and 83'. 

The various components of the propagator will be displayed with the LA and R contri- 
butions added up. 

For the "first propagator" (two pairs of replica indices coinciding) we have: 



Uo + 2j x ) kdk k + 2^ + J X kdki 



rx 1 fXl ^ dk q ^ i _ r xi dk d r xi dl 8 

Tdi k > 1 

X1 f dk 8 dl 8 \ 

k 8k I 81 ^ + x ^ ■ [ } 



X 



For the "second propagator" (one pair of replica indices coinciding) we have two different 
formulae according to the relative values of the variables: 



I rt 1 f xi \ dk 8 x x 1 f x i dk 8 r* dl 8 x , 

' 1 ' t ~ Wo + 2J x + 4Jt )Tdk Fk ' + 4 F * + L TdkhTW^ 



Finally, for the "third propagator" (all replica indices different) we find four different 
expressions, depending again on the order of the variables: 

1 fv 1 ft\ dk 8 xv , s 

' +7 / , s<y<f<l, (H3) 



QX,y 


— u x,t 




G x ' v 


— ^x,t 


— QV, X — 


G x ' v 


_ q x <v 


— QV, X — 

n x,x 
^s,t — 



J ^ Jx 



2J X AJ y J k dk 
*\ dk 8 



k ()k F X k ' V , x<t<y< 1- (H4) 



dk 8 



■ kgk *r , » <»<»<i. (us) 



x 1 f s 1 /-*\ dk 8 



o 



+- / +- I ) + 



p x , x 

k dk 
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f s dk d [* dl d ^ x x . , 

+ J. Ten J. 7**! ■ **••'<»■ < 116 > 

In Eqs. 110-116 F*' y is the LA kernel given by Eqs. 100-102 and F%f is the R kernel 
defined in Eq. 86. 

The various propagator components are continuous functions of all overlap variables, 
except when either or both lower indices take on the value 1 corresponding to coinciding 
replica indices. 

8 The far infrared region 

The formulae derived in the previous section give a complete solution for the propagators 
of the Ising spin glass near the transition temperature and in zero field. Expanding them in 
the limit p 2 3> uq 2 one can easily recover the results, derived by different means in Sec. 6, 
for the near infrared region. To extract information concerning the far infrared (p 2 ~ uq 2 
or p 2 <C uq 2 ) is considerably harder. The main difficulty consists in the extreme richness of 
behaviour shown by the various propagator components. Indeed, depending on the relative 
magnitude of the overlap variables we have such a large number of different analytic forms 
that any attempt to display an exhaustive set of results would be quite illusory. What we 
are able to do is to merely give a sample of the results which will illustrate the types of 
behaviours one encounters in the extreme long wavelength limit. 

Let us start with the simplest propagators, with all overlap variables above the break- 
point (the "single- valley propagators"). With some effort one finds from the complete 
formulae in Sec. 7 that for p 2 <C uq 2 , i.e. far below the upper edge of the small mass band, 
these components are given by the following simple expressions: 

GiT « J2 . CUT) 

G£? - ' (118) 

G%\%\ « j 2 • (H9) 

The remarkable fact about Eqs. 117-119 is that they are the same as the p 2 — > limits 
of the formulae in Eq. 58 which were derived in the opposite limit, p 2 3> uq 2 . Although 
there are some (explicitly known but complicated) corrections around p 2 ~ uq 2 , the order 
of magnitude of the propagators is still ~ 1/p 2 , which means that the formulae given in 
Eq. 58 can be regarded as a good representation of the single-valley propagators in the 
entire momentum range. This was first pointed out in [30], but the p~ 2 like IR behaviour 
of the single valley propagators had been know already from [27] . 

We are not aware of any a priori reason why any propagator component should behave 
the same way both in the near and in the far infrared (and, indeed, no such continuity is 
observed in the other components), so we believe this coincidence carries a physical message. 
We shall return to this point in the next section. 
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From the continuity of the single-valley propagators one would tend to infer the same 
for the transverse and the longitudinal combinations given in Eq. 59. This is certainly true 
for the transverse component since, as shown already in [28], G± = 1/p 2 is an exact result 
within the Gaussian approximation. As for the longitudinal component, its behaviour in 
the far infrared is less clear. The point is that the leading terms cancel in the particular 
combination Gy and next to leading terms are not controlled reliably by the truncated model 
that we are using here. Nevertheless, it is safe to say that the longitudinal propagator is 
less singular than 1/p 2 for p —> 0. 

It is interesting to note that in some important respects (masslessness, the ratios 3:§:1 
between the three propagator components, weaker singularity in Gy) the qualitative be- 
haviour of the single-valley propagators of Eqs. 117—119 is rather similar to that of the 
correlation functions predicted by the droplet theory [4,5], despite the obvious conflict be- 
tween the underlying physical pictures. 

Let us now turn to the other propagator components. For the diagonal propagator G x ' x 
we find: 



yX,X U / 1 1 



w 2 \xp 3 2x 2 p 1 



0<x<xi, p^O , (120) 



where p is the rescaled momentum introduced in Eq. 95. Although G 1 ' 1 is evidently 
continuous at x = X\, Eq. 120 does not match Eq. 117. The apparent contradiction is 
resolved by noting that the limit p — > is nonuniform: near x = x±, G x,x depends on 
combinations like {x\ — x)jp and {x\ — x) 1 / 2 /p. The limit p — ► is nonuniform also around 
x = 0, so for a fixed value of x £ (0, 1) Eq. 120 will hold only if the wavenumber is sufficiently 
small to satisfy both pCi and p <C [x\ — x). Similar remarks apply in all the cases that 
follow. 

For the propagator component G x ' x , x < z, we find 



G Z = -2 b-o^l if °<*<^<*i . (121) 



u_ t J_ 

w 2 \xp 3 2x 2 p 2 

Again, in Eq. 121 [z — x) is understood to be large compared with p. 
The component G^, < x < y < x\ tends to a constant 

The divergence of this constant for x —>■ x\ — is a signal of the IR divergence of the limit 
x = y = xi, as given by Eq. 118. For x < y < x\ and p — > 0, G x ^ v x tends to another constant 
which is too complicated to be recorded here. 

The last item in this category is G x ' x . For < x < x\ we find 

G T X = ^{\-14^) > p^°- ( 123 ) 

L ' x w 2 \xp 3 2x 2 p 2 J 

Now we list a number of results for the "third kind" of propagator components with 
neither of their lower indices equal to 1. 
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Let us start with G%>%. For a fixed x € (0, x\) and p^Owe find: 



1 



7 1 



QX,X 

x ' x w 2 \xp 3 2x 2 p 2 



, < x < x\ , p — > 0. 



(124) 



Next we consider G^, < x < y. For 0<x<y<xiwe find a qualitatively new type 
of behaviour: 

7 1 



rtx,y 



^X'X 









\xp z 








yxp^ 



2 x 2 p 
7 1 



exp 
exp 



x\ — x 
p 

y-x 



P 



, 0<x<y = xi , p^O, (125) 
, 0<x<y<xi , p^O. (126) 



In the long wavelength limit, p — > 0, Eqs. 125 and 126 become (5-like distributions. 
Distribution-like components were first identified among the spin glass propagators by two 
of us [48], see also [49]. Independently, Ferrero and Parisi put them to use in their recent 
analysis of the IR divergences in spin glasses [24]. 

Returning to the list of propagators we find that the component G^'V, < x < y < x±, 
is equal to the constant in Eq. 122 for < x < y = xi, and 

8u 1 S(x/y) 



rtx,y 
x,y 



Furthermore, 
G 

G 



x,x _ 


u 




x,y ~ 


w 2 


yxp^ 


x,x 


u 


{— 


y,y ~ 


w 2 


\yp 2 


x,x 


u 




y,y ~ 


w 2 





w 2 y 2 {y 2 -x 2 ) 5(1) 
5 1 



2 x 2 p 2 
7 

~2 



y, p 

7 1 





f 2 Xl ~ y ] 


) exp 






P . 



, < x < y < xi , p^O. (127) 

< x <y<Xi , p^O, (128) 
, y<x = xi , p^O, (129) 



exp 



x-y 



P 



, y <x <xi , p ->■ 0. (130) 



The remaining components depend on three different overlap variables: 



/~1X,X 

^y,z 



(~,x,y 
^x,z 



r<x,y 
^x,z 



(~ix,y 

^x,z 



(~ix,y 

^z,z 



fix,y 
^z,z 



rix,y 

W 2,2 



1 



1 1 



w 2 \ xp 3 2 x 2 p 2 

4u 1 S(x/xi 
w 2 xi(x 2 — x 2 ) 5(1) 

lQu p S(xjz) 
w 2 z 2 (z 2 — x 2 ) 5(1) 

8u 1 S(x/z) 



, 0<x<y,z<xi , p — > 0, 



exp 



w 2 z 2 {z 2 -x 2 ) 5(1) 
7 



exp 



P 

y - z 



P 



(131) 
(132) 

, 0<x<z<y = x±, p — > 0,(133) 
, 0<x<z<y<x 1 , p -> 0, (134) 



, 0<x<y = z = x\ , p — > 0, 
X\ — z 





(A 




\zp z 








\zp z 








yzp 3 



7 
z 2 p 
7 1 



exp 



exp 



-2 



xi — z 



P 



, z < x = y = x\ , p — > 0, (135) 



xi + x — 2z 



2 z 2 j5 2 



exp 



y + x — 2z 



P 



, z < x < y = xi , p — > 0, (136) 
, z<x<y<xi , p^0, (137) 
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and finally G^f goes to different constants as p — »■ 0, for x < y = z = x\, x < y < z = x±, 
and x < y < z < Xi, respectively. 

The components with one or two upper indices equal to zero are special. If the two upper 
indices are different and the smaller is zero, then the propagator component in question 
vanishes indentically. If the two upper indices coincide and vanish, then the leading infrared 
behaviour in the diagonal component is given by 

G i i = T — ^4 ' P^O. 138) 

As long as the lower indices s, t are larger than zero, the off-diagonal components and 

are the same as Eq. 138; for t — > (or s and/or t — ► 0) they vanish again. 

The strong, p~ 4 like singularity in a correlation function related to Eq. 138 was first 
noticed by Sompolinsky and Zippelius [50]. Several results concerning the far infrared region 
were published by us [27—30], by Goltsev [51], and by Ferrero and Parisi [24]. 

The compilation of Eqs. 117—138 above gives more details than any of these papers, 
nevertheless it is far from exhaustive. In a propagator component depending on, say, two 
overlap variables we may have, e.g., that one of them is of the order of x±, while the other 
is much less than even x\ , thereby defining a new mass scale. Depending on the value of 
the momentum relative to these new mass regions we may have further details in the IR 
behaviour. With the exception of Eq. 138 the results listed in this section have been derived 
under the assumption that all overlap variables are of the same order of magnitude as x\ 
and the momentum is much smaller than x\ . 



9 The physical meaning of the propagators 

In addition to their role as building blocks of the interacting field theory, the Gaussian 
propagators have also a direct physical meaning which we would like to discuss now. This 
will also allow us to give a physical interpretation to some of the results we have found so 
far. 

The propagators are related, as in any field theory, to some correlation functions. These 
correlation functions reflect the underlying assumption we have adopted about the structure 
of phase space, namely that it splits into a large number of equilibrium states with a 
hierarchical organization. 

The simplest correlation function one can define is the overlap of the spin-spin correla- 
tion function < SiSj > in valley a with the same in valley b: 

C ah (r) = — < s i s i >a< SiSj > h = q& + — ^ e^CUp) (139) 

i P 

where the distance r = r-j — between the sites is kept fixed as we sum over the lattice 
points i, and q a b is the overlap between states a, b, Eq. 31. 

In principle, C a b{r) as defined in Eq. 139 could depend on the concrete realization of the 
random couplings Jjj and also on the two states a and b. Following the considerations of 
Parisi [38] and also those of Mezard and Virasoro [52] one is led to the conclusion, however, 
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that C a b(r) is self-averaging (i.e. independent of the sample Jij for a large system) and 
depends on a and b only through the overlap q a b- Moreover, the Fourier transform 

Cab(p) = -^X)exp[i(rj ~ r i)p] < s i s j >»< s i s j >b ~ N <llb S p,o 

ij 

can be calculated in the Gaussian approximation via the replica formalism and turns out 
to be nothing but the diagonal component of the propagator: 

C ab (p) = G?£(p) (140) 

where the relationship between a, b and x is given by the inverse of Parisi's order paramater 
function: x = x(q ab ). The derivation of Eq. 140 follows the same steps that led Parisi to 
dx 

the identification of — as the probability distribution of overlaps [38], and can therefore 

be omitted here. The essential ingredient is the assumption about the states having the 
property of clustering. 

Correlation functions involving three or four states can be treated similarly: 

C abc (r) = -j- < s i s j >a< Si > b < Sj > c = q ab q a c + T7 H e" ipr C afec (p) (141) 



i P 



and 



1 1 

Cabcdi?) = < Si >a< Si >b< S i >c< S i >d= QabQcd + e ' 1PTC abcd{p) (142) 

i P 

are also self-averaging and depend on the states abc, resp. abed only through the overlaps. 
The Fourier transforms are given by the off-diagonal components of the propagators as 

Cabdv) = G?;«(p) (143) 

with x = x(q ab ), y = x(q ac ), z = max{x(^ a ), x(q bc )}, and 

C abcd (p)=Gfy z2 (p) (144) 

with x = x(q ab ), y = x(q cd ), z\ = m&x{x(q ac ), x(q ad )}, z 2 = max{x(g &c ), x(q bd )}. 

Eqs. 140, 143 and 144, first written up in [30], give a direct physical meaning to the 
various propagator components. We see now that the complicated behaviour of the propa- 
gators reflects the structure of correlation overlaps inside a state and also between different 
states. 

In particular, we can see that the propagator components with all overlaps above x\ = 
a?((/max) correspond indeed to correlation functions inside a single state, as anticipated in 
Sec. 6. In the near infrared region we found for these propagators the simple forms in 
Eq. 58, and the transverse and longitudinal combinations given in Eq. 59. These clearly 
suggest the idea of a massless phase with a transverse correlation function falling off like 
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l/r d 2 in real space and an exponentially decaying longitudinal correlation function with a 
characteristic length 

! ' (145) 



s \Z2wqx \/2t 

This reveals the meaning of the large mass: it is the inverse coherence length of the longi- 
tudinal fluctuations in a single state. 

In order for this interpretation to be consistent, the single-state propagators should not 
develop stronger infrared singularities in the far infrared region than the behaviour we have 
just seen in the near infrared. Recalling Eqs. 117, 118 and 119, valid for momenta far below 
the upper edge of the small mass band, we see that this is precisely what happens: the 
leading terms 3/p 2 , 3/2p 2 and 1/p 2 one finds for the single state propagators in the far 
infrared coincide with the p — > limit of the formulae (58), derived in the near infrared. It 
seems therefore that the long-distance behaviour of the single-state propagators is relatively 
mild: although the phase is massless, the infrared behaviour does not appear to be more 
violent than in the Heisenberg model, for example. In addition, as shown in [42], the 
result 1/p 2 for the transverse single-state propagator is, in fact, exact within the Gaussian 
approximation; it holds not only in the truncated model or near T c , but throughout the 
ordered phase. 

For a complete characterization of the correlations in the system we have to study 
also the interstate overlaps of correlation functions, with some or all replica overlaps going 
below x\. The behaviour one finds in this case depends radically on the momentum range 
considered. For momenta in the near infrared the correlation functions display a reasonable 
degree of complexity: p scales with the large mass, and the overall scaling power is 1/p 2 , as 
shown by Eqs. 60—69. This means that the qualitative features of the correlation overlaps 
between two states with a given overlap q < g max are similar to those inside a single state, 
except on the extremely long wavelength scales. The characteristic length beyond which 
qualitative changes occur is £' ~ 1/r, corresponding to the upper edge of the small mass 
band. This gives a meaning to the small mass scale. 

As we see from the results in the previous section, for wavenumbers in the range of the 
small mass or below, the correlation overlaps show very strong infrared divergences, going 
up to 1/p 4 . This violent IR behaviour poses a formidable challenge to spin glass theory. 



10 The first-loop corrections above 8 dimensions 

In this section we apply the previous results for the calculation of the first short-range 
corrections to the equation of state and to the mass spectrum. The treatment follows our 
previous works [15,16] closely. For the time being we assume that the dimensionality is 
sufficiently high for the ordinary loop-expansion to work: as we shall see shortly, this means 
d > 8. 

Our starting point is Eq. 18 where we have to remember that q a p is, in principle, the 
exact order parameter, but r is the reduced temperature measured relative to the mean field 
transition temperature. Let r c be the critical value of r where the exact q a p vanishes. Now 
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we can write up Eq. 18 at r c with the help of Eqs. 40 and 41 which are valid for q a g — ► 0, 
and subtract the resulting equation from (18). We get: 



2u 

2tq aB + w(q 2 ) a3 + —qls + (146) 



1 1 



+ 'z (27 



( d d p W ( G °n,-lP (P) ~ + 2u 9a/3 ( Ga/?,a/?(p) - — 
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where £ = r — r c is the "true" reduced temperature. The critical value of r works out to be 

r c = - z (w 2 h ~ ul 2 ) , (147) 

where 

f d d v 1 

J '= / TT^d- ' Z = 2,4,6,... . (148) 
V|p|<i {2-K) d p l 



(For / > d an IR cutoff at p > V2£ is understood in the integrals //.) 

Since r c is of 0(1/ z) we can replace r by £ in the propagators appearing in the loop 
terms. These propagators satisfy the set of equations (36) — (39), with r replaced everywhere 
by £. In the dimensionality range considered here the dominant contributions to the loop 
integrals come from the neighbourhood of the upper cutoff. In this momentum range the 
"large-p expansion" results of Eqs. 42—44 are valid. As the mean field part of the equation 
of state goes up to 0(q 3 ), we have to stop at the same order also in the loop terms. This 
means we have to go to 0(q 5 ) in G a7j/ g 7 in Eq. 43, but have to stop at 0(q 2 ) in G a p^ a p, 
Eq. 42, which is already multiplied by q a p. Substituting these approximate propagators into 
Eq. 146, performing the summation over the replica index 7 and collecting the coefficients 
of q a p, (q 2 )a/3, and q^g, we end up with an equation which is of the same form as the 
mean field equation of state with £, w and u replaced by some new, renormalized coupling 
constants £, w and u. The mapping works out as follows: 

t = £^l-i[4ti; 2 / 6 + 12ti; 2 £/ 8 -2n/4-4n£/ 6 ]^ 

w = w^l-^-[2w 2 lQ + 12w 2 tIs]j (149) 

(1 12-u; 4 \ 
l--[12w 2 I 6 — I 8 -6uI 4 ]\ ■ 

In addition to the above, some new type of couplings are also generated by the loop 
terms. These would be corrections to those quartic couplings which we discarded already 
at the very beginning. Since they enter only subleading terms in all quantities which we 
are interested in, we can safely ignore them. 

With this we have mapped back our loop-corrected theory onto MFT. The order param- 
eter will be of the same form as Parisi's order parameter function given in Eqs. 27 and 28, 
with r replaced by £ and w, u replaced by w, u, respectively. A more detailed discussion of 
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the effects of these replacements can be found in [15,16], here we merely note that the slope 
w/2u of q(x) decreases, while the breakpoint x\ increases with decreasing dimensionality. 
This means that replica symmetry breaking effects are enhanced by the loop corrections, in 
complete agreement with the conclusions drawn by Georges, Mezard and Yedidia from the 
1/d expansion [14]. 

We can arrive at the same conclusion by considering the excitation spectrum of the 
system. In order to be able to calculate the renormalized masses we need first the corrections 
to the self-energy matrix. At the first loop level these are typically bilinear combinations of 
the various propagator components, integrated over the loop momentum. For our present 
purposes we may disregard "wave function renormalization effects" and evaluate the self- 
energies for zero external momentum. For d > 8 the propagators entering the calculation 
can again be replaced by their "large-p" approximations from Eqs. 42—44 where we can now 
stop at 0{q 2 ). Performing all the necessary replica summations and momentum integrations 
we find at the end that the self-energies to 0(l/z) are precisely of the same form as at the 
zero loop level, Eq. 85, with the replacement of r, w and u by the renormalized coupling 
t, w and u, respectively. This shows the full consistency of our scheme and guarantees 
that important qualitative features of the spectrum (e.g. gaplessness) are not violated. The 
renormalized masses themselves can then be obtained from their mean-field counterparts, 
with the replacement of the bare couplings by the tilded ones. Considering Eq. 109 we 
can then see that the centre of the large mass band will change very little but, due to the 
renormalization of the u coupling, the bandwidth will grow considerably as the dimension 
d decreases. The same is true for the upper edge of the small mass band. These are further 
signals of growing RSB effects (without RSB the spectrum would consist of two points). 

We can sum up these findings in the following way: Parisi's RSB field theory has proved 
to be perturbatively stable above d = 8 dimensions. The loop corrections do not change 
the qualitative features of the theory, they just shift the constants by a small amount, and, 
at least to the low order regarded here, they even enhance RSB effects. 

As we approach d = 8 from above, however, the mapping set up in Eq. 149 becomes 
singular: the integral Is (cut off at the lower end at p = y/2t) blows up. This causes the most 
serious problem in u where Ig appears alone, without any additional t factor multiplying it. 
As it stands, u simply does not make sense below d = 8: for t — ► it diverges like t( rf ~ 8 )/ 2 . 
The next section is devoted to the resolution of this paradox. 

11 Between 6 and 8 dimensions 

The formal reason for the mapping (149) becoming singular is easy to understand: whereas 
for d > 8 the loop integral Is is dominated by the contribution of momenta around the UV 
cutoff, for d < 8 the main contribution comes from the lower end where the momentum is 
of the same order of magnitude as the large mass. Then the large-p expansion is not a good 
representation of the propagators any more. In particular, as the p~ 8 singularity comes 
from the expansion of G a7;7j g, we have to treat this component more carefully. (The other 
propagator, G a p^ a p, in the u-loop term in (146) does not cause problems for d > 6, so we 
can continue to use the large-p expansion result Eq. 42 for it.) 
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Working out the replica summation in Eq. 146 we get in the continuum limit: 

E G ^nP = ~ T dyG™ - xG x £ - 2 f dyG^J , an(3 = x . (150) 

Considering that x < X\ ~ t <C 1, it is clear that the largest term in Eq. 150 comes from 
the upper end of the last integral, so we can again use the "innermost block approximation" 
here: 

E G <*7,7/3 ~ ~ 2G x,T ■ ( 151 ) 

Now one can check that for d > 6 the contribution of the far infrared region to the 
integral of Eq. 151 over the momentum is still negligible, so we can use Eq. 62, valid in the 
near infrared, for the propagator in Eq. 151. Then substituting Eqs. 62 (taken at y = x±) 
and 42 into Eq. 146 and collecting the coefficients of q(x) and q 3 (x), respectively, we arrive 
finally at a new set of renormalized coupling constants which are, in fact, very similar to 
those in Eq. 149. The only difference is that u is now given by 
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U = U+ 12 ^-Tl—TT^ \~o , (152) 

z J {2ir) d p 4 (p 2 + 2t) 2 v ' 

where 2wq\ in the denominator of the integrand has been replaced by its zeroth order value 
2t. In the dimension range 6 < d < 8 regarded here the integral in Eq. 152 is well-defined 
both for small and for large momenta, so to leading order in t we can send the upper cutoff 
to infinity. The effective coupling u then becomes 

t (d-8)/2 

u = u + const • . (153) 

for fixed t and z very large (i.e. for temperatures not very close to T c and for very long range 
forces) the shift in the quartic coupling is small, and we are back to the previous situation: 
the loop-corrections do not significantly alter the mean field results. In the opposite limit, 
however, where z, albeit large, is kept fixed and the temperature is allowed to go arbitrarily 
close to T c , the bare u becomes negligible compared with what was supposed to be a small 
correction. It is clear that under these circumstances the ordinary loop-expansion breaks 
down. 

Standard power counting arguments tell us that the upper critical dimension of the 
model defined in Eq. 2 is d c = 6. It would therefore be perfectly normal to find an in- 
frared breakdown of the perturbation expansion in six dimensions. Why we should have IR 
problems already in d = 8, however, demands explanation. 

To understand the origin of the problem, we have to realize that the term that blows up 
in 8 dimension is a one-loop contribution to the 4-point function at zero external momenta. 
It is evident that this quantity (the "box graph" ) should be singular at the critical point in 
d = 8. 

It should be emphasized that this effect is not a spin glass peculiarity. Many-point 
functions at exceptional momenta (where power counting arguments fail) can blow up at 
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T c in high dimensions in any theory. For example, the one-loop correction to the sixth 
derivative of the thermodynamic potential with respect to the magnetization in an ordinary 
ip 4 theory (the triangle graph) will be singular at T c in d = 6, i.e. above the upper critical 
dimension of that model. This six-point function is related to a higher order nonlinear 
susceptibility, which is normally of little interest, so the problem is rarely mentioned. The 
triangle graph (and all the higher polygons that have similar IR singularities in higher 
and higher even dimensions) enter, however, as insertions also in some of the high order 
loop-corrections to physical quantities such as the two-point function (i.e. the inverse linear 
susceptibility) which are regularly discussed and claimed to be free of IR problems above the 
upper critical dimension. One may wonder whether these IR-singular many-point insertions 
do not make also the 2-point function singular. They do not: the singularity in the many- 
point functions occurs only at exceptional (zero external) momenta, so when these graphs 
appear as insertions in higher order diagrams and finite momenta flow through them their 
singularity will be suppressed by the loop integrals. Therefore, the isolated IR singularities 
appearing in the many-point functions do not proliferate in higher orders and do not destroy 
the loop-expansion for the quantities which one is usually interested in. 

The situation in spin glass field theory is quite similar to the above, except that here 
the many-point functions do appear directly in the coefficients of the Taylor expansion of 
the order parameter in the variable x. In particular, the 4-point function appears in the x 3 
term in the equation of state, hence in the slope of q(x), so if we want to calculate q(x) to 
linear order we have to learn how to handle this difficulty. 

Since above 6 dimensions the singularity will not proliferate in higher orders, it is evident 
what we have to do. We have to absorb this isolated singular term into the mean field part 
of our equation and treat the rest as perturbation. This way we can save the loop-expansion 
between 6 and 8 dimensions in a reorganized form. 

The starting point in this new expansion is an effective mean field theory which is 
exactly of the same structure as what we have discussed so far but with the bare coupling u 
replaced by u everywhere. Since the remaining loop-corrections do not qualitatively modify 
the results, in the following we restrict ourselves to the discussion of this effective MFT. 
Moreover, we will be interested only in the fluctuation dominated regime where 

*(d-8)/2 

> 1 , (154) 

so the bare u can even be completely omitted. 

The nontrivial renormalization of the four-point coupling has a profound consequence 
for the structure of the theory. To fully appreciate its significance we have to go back for 
a moment to the original MFT. As first stressed by Fisher and Sompolinsky [53], scaling is 
badly violated in Parisi's MFT: not only the so called hyperscaling laws break down, but 
also those that do not explicitly contain the dimension d, such as e.g. (38 = 1 — a/2 + 7/2. 
As a matter of fact, it is not even possible to unambiguously assign a critical exponent to 
some physical quantities. For example, the maximum of the order parameter function scales 
as qi ~ t, but Sompolinsky's susceptibility anomaly [54] A = ^ (qi — Jq 1 dxq(x)j which is 
a particularly useful measure of ordering in the presence of an external magnetic field goes 
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as A ~ t 2 . The existence of two "mass scales", i.e. two characteristic lengths, one diverging 
as i^ 1 / 2 the other as t -1 , is another example of this ambiguity. Now the most important 
effect of the replacement u — ► u below 8 dimensions is that upon approaching the upper 
critical dimension d = 6 from above scaling gets gradually restored [53], [15,16]. 

Indeed, consider e.g. the slope c = w/2u of q(x). With u ~ t d / 2 ~ 4 /z, valid in the 
fluctuation dominated regime, the slope becomes c ~ zt 4 ~ d l 2 . For d — > 6 + c is therefore 
~ t. On the other hand, the plateau of q(x) is independent of u, so it remains q± ~ t. As for 
the breakpoint x\ = 2ut/w 2 , it becomes ~ t d l 2 ~ 3 jz. Note that for d —> 6 + the breakpoint 
becomes independent of the temperature (but remains small, of the order of 1/z). What all 
this amounts to is that upon approaching d = 6 from above the order parameter function 
becomes of the form 

q(x) = tPf(x) , (155) 

where ft = 1, and f(x) is independent of temperature. The form (155) is necessary for all 
the meaningful combinations one can form from q(x) to scale with the unique exponent (5. 

Similar remarks apply for the mass scales. The center of the band of the large eigenvalues 
remains 2t, but as shown by Eq. 109 the band width 4t 2 u/3w 2 varies with d as ~ 

t d/2-2 /z 

and becomes ~ t/z as d —> 6 + . The same is true for the upper edge 2ut 2 /w 2 of the band of 
small eigenvalues. This means that instead of two mass scales having different temperature 
exponents, we have a single exponent v = 1/2 for t — ► 6 + . Note however that the separation 
of scales still remains in the form of a numerical difference: the small eigenvalues and the 
band width of the large eigenvalues is down by a factor 1/z compared with the center of 
the large eigenvalues. 

To make contact with the work of Green, Moore and Bray [55] (who were, in fact, the 
first to note the role of the renormalized four-point coupling in a particular instance), we 
consider now the effect of an external field h on the order parameter q{x). Although we 
have avoided the problem of the field in this paper, it is easy to show that the replacement 
u — > u will work also in the presence of the field. Therefore the well-known results for the 
field dependence of q(x) can be readily taken over from ordinary MFT, and one can see, 
in particular, that the AT line [31] will be given by h 2 AT = 4St 3 /3w 3 . The standard mean 
field result h 2 AT ~ t 3 will then be replaced by h 2 AT ~ t d l 2 ~ x jz for 6 < d < 8, becoming 
h 2 AT ~ t 2 /z as d — ► 6 + . This coincides with the result found by Green et al. [55] which they 
obtained from the zero of the replicon self-energy in a one-loop calculation in the disordered 
phase. The agreement between these two independent calculations is a testimony for the 
consistency of the effective MFT. Green et al. also suggested that the shifted AT exponent 
might be exact. This must indeed be true for all the exponents predicted by the effective 
MFT, because, as argued above, the singularity found in the four point function will not 
proliferate and the higher loop corrections will not be more singular then the one-loop 
contribution worked out here. 

When we were talking about the limit d — > 6 + we always meant to remain slightly 
above d = 6. In exactly six dimensions the other corrections, disregarded in the effective 
MFT also become singular and produce the usual logarithms appearing at the upper critical 
dimension. The little we know about the range d < 6 is the subject of the next section. 
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12 The first step below 6 dimensions 



As we cross d = 6 infrared divergences start to plague the whole perturbation expansion. 
Now they do proliferate: every order will be more singular then the previous one. In this 
situation one normally turns to renormalization group (RG) methods to reorganize the 
perturbation expansion. The expansion in the distance e from the upper critical dimension 
proved particularly successful in the context of ordinary critical phenomena. The application 
of field theoretic RG methods to spin glasses was pioneered by Harris, Lubensky and Chen 
[32]. Building on [33], Green calculated spin glass critical exponents to 0(e 3 ) in [34]. The 
structure of the RG in spin glasses is poorly understood, however. No systematic RG 
analysis has ever been made in the condensed phase, and the negative result of Bray and 
Roberts [56] who failed to locate a stable fixed point in a field below d = 6 remains an 
enigma to resolve. 

Neither do we have a magic key to the riddle. What we are able to do here is to take just 
the very first step toward an e-expansion in the ordered phase of spin glasses, without the 
solid basis of an RG analysis. The findings in the previous section provide some important 
clues. 

The first of these is that the bare quartic coupling u can be dropped altogether — the 
cubic theory builds up its own quartic coupling at the first loop level. This was already 
observed by Bray and Moore [57] who, working around a replica symmetric starting point, 
run into an AT-like instability, however. Secondly, with the restoration of scaling the 
breakpoint x\ of the order parameter becomes small, of 0(1/ z). If there exists an e- 
expansion at all, this should translate into x\ = 0(e) below d = 6. But if x\ is small, so 
is the whole region where RSB occurs, hence to leading order we can apply the innermost 
block approximation in all replica summation. Thirdly, the upper edge of the small mass 
band must also become 0(e) , so to leading order we can always use the near infrared 
propagators. 

With all these simplifications the calculation of some of the critical exponents to 0(e) 
becomes quite feasible, for the exponent [3 which we present here as an illustration, it 
becomes almost trivial. 

In order to obtain (3, it is sufficient to calculate q\. Let us regard the equation of state 
as a polynomial in x. Then a little reflection shows that the maximum of q(x) is determined 
by the linear terms in x. (The slope at x = would be determined by the x 3 terms.) 

Where do we get linear terms in Eq. 146? The first term is such. In the second term 
we have 

(q 2 )*p = - f X dyq 2 (y)-xq 2 (x)-2q(x) [' dyq(y) , a D (3 = x , (156) 

JO Jx 

which provides a linear term coming from the last integral 

(q 2 U^-(2- Xl )q iq (x) + 0(x 3 ) . (157) 

With the bare u omitted the only term we are left with is the w-loop. Using the innermost 
block approximation (151) for the replica summation, and selecting the linear term from 
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the near infrared approximation (62) for the propagator appearing here we have 

Ek^-^=-l!?^ + ^ 3 ) , (158) 
4^3 V ^ P 4 J p 2 (p 2 + 2t)i 1 ; y J 

where, in view of this loop term being of 0(e), we used the zeroth order result wq\ = t. 

Collecting Eqs. 156—158 and evaluating the integral of Eq. 158 in d = 6 for t — ► 0, we 
find that the equation of state to first order in x is: 



^2t - 2wq 1 + wxiqx - Aw ^ m y/2tj q ( x ) + Q(x 3 ) = , 



(159) 



where K 6 = 

Writing wqi = At 13 = (1 + A^> -\ )t(l + [3^ hit H ) where A^ and f3^ are the 

0(e) corrections to the amplitude A and to the critical exponent (3, respectively, we finally 
obtain 

/3« = ^ . (160) 

From [33,34] we know that at the fixed point of the cubic coupling the combination in 
Eq. 160 is e/2. Thus we get 

/3 = l + | + 0(e 2 ) 

which is in accord with the exponents calculated at T c [34]. 

If the interpretation of the IR logarithm coming from the loop-integral as a correction 
to is correct then at the next order the coefficient of In 2 t must be half of the coefficient 
of the first order log. We have evaluated the two-loop correction to the equation of state to 
In 2 1 order and found that the above condition of exponentiation is fulfilled, indeed, provided 
(160) is chosen e/2, its fixed point value. 

Although this provides an important check on our first order calculation and raises 
the hope that the e-expansion can be consistently carried through also below T c , we have 
to point out that at the second order we find not only the ln 2 i, hit terms necessary for 
scaling, but also infrared divergences like p~ 6 . In higher orders we will clearly find even 
stronger infrared powers, due to the presence of zero modes and soft modes in the system. 
Though we do notice some cancellations, a systematic method (similar to the exploitation 
of rotational symmetry through the use of Ward identities in the 0(n) model, or the use of 
the small masses as an infrared cutoff) is still to be found in spin glass theory. Work is in 
progress in that direction. 



13 Summary 

Let us briefly recapitulate the main points in the paper. Having set up the field theoretic 
formalism, we devoted a long discussion to the free, quadratic fluctuations about an equi- 
librium solution with Parisi's ultrametric symmetry. We showed that for momenta higher 
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than any characteristic mass, i.e. for wavelengths shorter than any characteristic distance 
in the system, the propagators describing these fluctuations can be obtained even without 
having to specify the concrete symmetry breaking pattern. 

In the long but not extremely long wavelength limit, that we called the near infrared, we 
obtained simple approximate forms (rational fractions with poles at zero and at the large 
mass) for the propagators. We also found that in this region the propagators show a simple 
and explicit functional dependence on the order parameter: they are polynomials, going up 
to the quartic order, in q(x). 

We also displayed the group theoretic and algebraic techniques one needs to block- di- 
agonalize any ultrametric matrix. The additional simplification occuring in the kernel of 
these block- diagonal forms in the case of the Gaussian propagators near T c then allowed 
us to obtain exact, closed expressions from which the far infrared behaviour could be de- 
termined. The dependence of the propagators on q(x) is much more complicated in this 
region than in the near infrared, in particular, the high IR powers we found for extremely 
small wavenumbers are intimately related to q(x) being a linear function for small x. 

The physical meaning of these results could be understood on the basis of the relation- 
ships we established between the propagators and some intra- and intervalley overlaps of 
spin-spin correlation functions. In particular, we learned that transverse fluctuations are 
long-ranged also inside a single phase, while the longitudinal fluctuations have a finite co- 
herence length, given by the inverse large mass. Overlaps of correlations between different 
phase space valleys were found to start to qualitatively deviate from those inside a single 
phase around a distance given by the inverse small mass. Most remarkably, inside a single 
phase the extreme long wavelength behaviour was found to be a smooth continuation of the 
behaviour on the intermediate, near infrared scales, whereas the overlap between correlation 
functions in two distant valleys shows a markedly different, more singular behaviour in the 
far infrared than in the near infrared. 

Having understood the structure of correlations in the various regimes we proceeded to 
apply our propagators in the calculation of the first loop corrections. We showed that for 
d > 8 the theory maps back onto MFT with small, numerical shifts in the coupling constants, 
which demonstrated the stability of Parisi's MFT against short range corrections, at least 
to the low order investigated here. 

For 6 < d < 8 we had to reorganize the loop expansion in order to incorporate a divergent 
correction to the quartic coupling. This led us to an effective MFT with exactly calculable 
but dimension dependent critical exponents, and enabled us to follow the development of 
the theory with decreasing dimensionality towards a form reached at d = 6 + , with scaling 
restored but a nontrivial replica symmetry breaking pattern preserved. 

Entering the range d < 6 we noted that both the breakpoint x\ and the small masses 
become of the order of e . This simplified the calculation tremendously, and allowed us 
to readily identify the term that yields the 0{e) correction to the critical exponent of the 
order parameter, and even to check exponentiation at the next order. 

We concluded the paper by pointing out the challenge posed by the appearence of 
infrared singularities in the expansion and by hinting at the methods by which this challenge 
may be met in the future. 
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